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Abstract. The set of all group relaxations of an integer program 
contains certain special members called Gomory relaxations. A 
family of integer programs with a fixed coefficient matrix and cost 
vector but varying right hand sides is a Gomory family if every 
program in the family can be solved by one of its Gomory relax- 
ations. In this paper, we characterize Gomory families. Every TDI 
system gives a Gomory family, and we construct Gomory families 
from matrices whose columns form a Hilbert basis for the cone 
they generate. The existence of Gomory families is related to the 
Hilbert covering problems that arose from the conjectures of Sebo. 
Connections to commutative algebra are outlined at the end. 



1. Introduction 

Given an integer dxn matrix A and a cost vector c G Z n , we consider 
the family IPa, c of all feasible integer programs of the form 

IPA,c{b) '■= min {c ■ x : Ax = b, x G N n } 

as the right hand side vector b varies. The matrix A is assumed to have 
rank d and cone(A), the cone generated by the columns of A, is assumed 
to be pointed. We also assume that {x G IR> : Ax = 0} = {0} which 
guarantees that all programs in IPa, c are bounded. 

In 0, Gomory defined the group relaxation of IPA, c {b): 

min {da ■ Xa : A a x a + A^x s — b, x„ > 0, (x a , x & ) G Z n } 

where A G is the optimal basis of the linear relaxation of IPA, c {b) and 
non-negativity restrictions on the optimal basic variables x a have been 
dropped. The cost vector c d is the restriction of (c — c^A^A) to the 
components indexed by the complement of a. The extended group relax- 
ations of IPA, c (b) introduced by Wolsey are the 2^ relaxations obtained 
by dropping non-negativity restrictions on each subset of the variables 
in x (j The set of (/-dimensional simplicial cones cone(A a ), as A a 



varies over the optimal bases of the LP-relaxations of all programs in 
IPa,c, triangulates cone(A). This is called the regular triangulation of 
cone(A) with respect to c and is denoted as A c . The collection of all 
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sets a indexing the full dimensional cones of A c along with all their 
subsets form a simplicial complex on {1, . . . , n} which we also call A c . 
In this paper we consider the set of all group relaxations of IPa, c (P) 
obtained by dropping non-negativity restrictions on the variables in- 
dexed by each face of A c (Definition |2.6| ). This is a larger set of group 
relaxations for IP^ c (b) than the set of extended group relaxations of 
Wolsey. We show that these are precisely all the bounded group re- 
laxations of IPA, c {b) (Theorem |2.7| ). Among these group relaxations, 
the easiest to solve are those indexed by the maximal faces of A c . We 
call these the Gomory relaxations of IPa, c {o)- The family of integer 
programs IPa, c is called a Gomory family if all its members can be 
solved by one of their Gomory relaxations. 

Theorem [2.171 characterizes Gomory families. This theorem is a con- 
sequence of re-casting algebraic results on toric initial ideals in terms 
of group relaxations of integer programs |^|, [14|. These algebraic re- 



sults come from Grobner bases methods in integer programming fl9 |. 
No familiarity with these techniques is assumed in this paper. In Sec- 
tion 3 we relate Gomory families to total dual integrality (TDI-ness). 
Theorem |3.3| shows that yA < c is a TDI system if and only if the reg- 
ular triangulation A c is unimodular. This leads to Corollary |3.4| that 
if yA < c is TDI then IPa,c is a Gomory family. 

In Sections 4 and 5 we exhibit general classes of Gomory families. A 
matrix A is said to be normal if its columns form a Hilbert basis for 
cone(A). In Section 4 we introduce A-normal matrices which form a 
proper subset of normal matrices. Theorem [4.7] shows that every A- 
normal matrix A gives rise to a Gomory family IPa,c- While we do not 
know if every normal matrix gives rise to a Gomory family, we show that 
for small values of d every regular triangulation of cone(A) supports 
Gomory families (Theorem |5.5| ). Gomory families induce special covers 
of NA, the semigroup generated by the columns of A, which relates their 
existence to the Hilbert cover questions found in 0, and |T7| . 

Throughout this paper we consider triangulations of cone(A). We 
wish to point out that a one-dimensional face of any such triangulation 
must be generated by a column of A. All computations in this paper 
rely on the connections of this material to commutative algebra as 
described in Ullf , [[0| and [^l]]. The relevant connections and 



codes are described briefly in Section 6. 
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2. GOMORY FAMILIES 

In this paper, we fix a matrix A G Z dxn of rank d and a cost vector 
c6Z" and consider the family IPa, c of all integer programs 

IPa,JJ>) ■= min {c ■ x : Ax = b, x G N n } 

as b varies in the semigroup NA := {Au : u G N n } C Z d . This 
semigroup is contained in the intersection of cone(A) := {Ax : x G 
M> }, and ZA := {^4z : z G Z n }, the lattice generated by the columns 
of A. We may assume without loss of generality that ZA = Z d . 
The feasible linear programs from A and c are of the form 

LPA, c {b) '■= min { c ■ x : Ax = b, x > 0} 

where b G cone(A). We denote this family as LPa, c - For a C {1, . . . ,n}, 
let y4 CT be the submatrix of A whose set of column indices is a. 

Definition 2.1. For a C {1, . . . ,n}, cone(Acr) is a face of the regular 
subdivision A c of cone(A) if and only if there exists a vector y G M d 
such that y ■ cij = Cj for all j G cr and y ■ Oj < Cj for all j o~. 

Remark 2.2. The regular subdivision A c is gotten by taking the cone 
in R d+1 generated by the lifted vectors (aj,Cj) G M d+1 where is the 
i-th column of A and q is the i-th component of c, and then projecting 
the lower facets of this lifted cone back onto cone(A). (See 0].) 

We assume that c is generic, which means that A c is a triangula- 
tion of cone(A). All cost vectors except those lying on a finite set of 
hyperplanes of M™ are generic 0. Using cr to label cone(A a ), the trian- 
gulation A c can be denoted as a set of subsets of {1, . . . , n}. This set 
of sets is closed under inclusion since A c is a simplicial complex, and 
hence it is specified completely by its maximal elements. For a vector 
x G R™, let supp(x) := {i : Xi ^ 0} denote the support of x. The signif- 
icance of regular triangulations for linear programming is summarized 
in the following proposition. 

Proposition 2.3. |2(| Lemma 1.4] An optimal solution of LPA, c {b) is 
any feasible solution x* such that supp(x*) = r where r is the smallest 
face of the regular triangulation A c such that b G cone(A T ). □ 

Proposition [O] implies that a C {1, . . . , n} is a maximal face of A c 
if and only if A a is an optimal basis for all LPA, c {b) with b in cone(A a ). 
Given a polyhedron P C M n and a face F of P, the normal cone of F at 
P is the cone N P (F) := {to G M™ : u-x 1 ><J'X, W G F and x G P}. 
The set of all normal cones of P form the normal fan of P in IR™. 
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Proposition 2.4. The regular triangulation A c of cone(A) is the nor- 
mal fan of the polyhedron P c := {y G M d : yA < c}. 

Proof. The polyhedron P c is the feasible region of max {y ■ b : yA < 
c, y G the dual program to LPa iC {o)- The normal fan of P c is 
supported on corie(A), i.e. the union of the normal cones of P c is 
cone(A), since this is the polar cone of the recession cone {y G M. d : 
yA < 0} of P c . Suppose b is any vector in the interior of a maximal 
face cone(A a ) of A c . Then by Proposition pT3| , LPa, c (1>) has an optimal 
solution x* with support a. The optimal solution y to the dual of 
LPA,c{b) satisfies y ■ a,j = Cj for all j G o and t/ • a j < Cj otherwise, 
by complementary slackness. Since a is a maximal face of A c , in fact, 
y ■ aj < Cj for all j G' a. This shows that y is unique, and cone(A a ) 
is contained in the normal cone of P c at the vertex y. If b lies in the 
interior of another maximal face cone(A T ) then y', the dual optimal 
solution to LPA <c {.b) satisfies y' ■ A T = c T and y' ■ A f < c f where r ^ a. 
Hence y' is distinct from y and each maximal cone in A c lies in a distinct 
maximal cone in the normal fan of P c . Since A c and the normal fan of 
P c have the same support, they must therefore coincide. □ 



Corollary 2.5. The polyhedron P c is simple if and only if the regular 
subdivision A c is a triangulation of cone(A). □ 

Regular subdivisions were introduced in || and have since been stud- 
ied from various points of view. They play a central role in the algebraic 
study of integer programming (|19|], [pD|). We use them here to define 



group relaxations of IPa, c {P). 

A subset t of {1, . . . , n} partitions x = (xi, . . . , x n ) as x T and x f 
where x T consists of the variables indexed by r, and x f the variables 
indexed by the complementary set f . Similarly, the matrix A is par- 
titioned as A = [A T ,Af] and the cost vector clS C — f C T j Cf ). If a is a 
maximal face of A c then A a is nonsingular and Ax = b can be written 
as x a = A^^b — A & x & ). Then c ■ x = c CT (v4~ 1 (fe — A & x & )) + c„Xa = 
c a A~ x b + (eg. — c a A~ 1 Ag)xg. Let c„ := c„ — CaA^A^ and for any face 
r of a, let c f be the extension of c s to a vector in M' T ' by adding zeros. 

Definition 2.6. The group relaxation of the integer program JP^ )C (6) 
with respect to the face r of A c , denoted as G T (b), is the program 

min {c f ■ x f : A T x T + A f x f = b, x f > 0, (x T , x f ) G Z n }. 

The group relaxation G T (6) solves IPA, c (b) if its optimal solution is 
non-negative. These relaxations contain among them the usual group 
relaxations of IPA, c {b) found in the literature. The program G a (b) 
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where A a is an optimal basis of the linear relaxation LP^db) is pre- 
cisely Gomory's group relaxation of IPA, c (b) 0. The set of relax- 
ations G T {b) as r varies among the subsets of this a are the extended 
group relaxations of IPa,c(P) defined by Wolsey p4|. Since G A c , 
= IPA, c (b) is a group relaxation of IP^db) by Definition 



and hence JP J 4 )C (6) will certainly be solved by one of its extended group 
relaxations. However, it is easy to construct examples where a group 
relaxation G T (b) solves IPA, c {b), but G T (b) is neither Gomory's group 
relaxation of IPA, c (b) nor one of its nontrivial extended Wolsey relax- 
ations (see Example [2.16| ). Theorem |2.7| will show that the relaxations 



in Definition |2.6| are precisely all the bounded group relaxations of 
7P J 4 ;C (6). Hence this definition considers all the group relaxations of 
each integer program in IPa, c that can possibly solve the program. 

For our purposes it is convenient to reformulate G T (b) as follows. Let 
B G ^ nx ( n - d ) be any matrix such that the columns of B generate the 
(n — (i)-dimensional lattice C = {x G Z n : Ax = 0} C Z" and let u be 
a feasible solution of IPA, c {b)- Then 



IPA,c{b) = min {c ■ x 
= min {c ■ x 
= min {c ■ x 



Ax = b, x G N n } 
x = u (mod C), x > 0} 
x = u- Bz, x>0, z G Z n - d } 
The last problem is equivalent to min {c ■ (u — Bz) : Bz < u, z G 
Z n ~ d } and hence, IPA,c{b) is equivalent to the problem 

(1) min {(-cB)z : Bz < u, z G Z"~ d }. 

There is a bijection between the set of feasible solutions of ([]]) and the 
set of feasible solutions of IPA, c {b) via the isomorphism z i— ► u — Bz. 
In particular, G W n ~ d is feasible for ([[]) since it is the pre-image of u. 

Let 7r T be the projection map from IR n — > R' T ' that kills all coordi- 
nates indexed by r. If B T denotes the |f| x [n — d) submatrix of B 
obtained by deleting the rows indexed by r, then we denote by C T the 
lattice tt t (C) = {B f z : z G Z n ~ d }. It can be deduced from ^ that 
the group relaxation G T (b) is equivalent to the lattice program 

min {c f ■ x f : x f = tc t (u) (modC T ), x f > 0} 

which can be reformulated as above to be 

min {(-c f B f ) ■ z : B f z < tt t (m), zeZ^}. 

Since c f = ir T (c — c^A^A) for any maximal face a of A c containing 
r, and the support of c — c a A~ l A is contained in f, we get that c f B T = 
(c — c u A~ 1 A)B = cB since AB = 0. Hence G T (6) is equivalent to 

(2) min{(-cB)-z : B f z < n T (u), zeZ n ~ d }. 
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The feasible solutions to ([!]) are the lattice points in the rational poly- 
hedron P u := {z G W 1 ^ : Bz < u} and those to (@) are the lattice 
points in the relaxation P^ := {z G W l ~ d : B T z < n T (u)} of P u ob- 
tained by deleting the inequalities indexed by r. In theory, one could 
define group relaxations of IPa, c (P) with respect to any r C {1, . . . , n}. 
The following result justifies Definition 



Theorem 2.7. The group relaxation G T (b) o//Pa jC (6) has a finite op- 
timal solution if and only if r C {1, . . . , n} is a face of A c . 

Proof. Since all data are integral it suffices to prove that the linear 
relaxation min {(—cB)z : z G P^} is bounded if and only if r G A c . 

If r is a face of A c then there exists y EM. d such that yA T = c T and 
yA f < Cf. Using the fact that A T B T + A f B T = we see that cB = 
c T B T + c f B f = yA T B T + c f B f = y(-A f B f ) + c f B f = (c f - yA f )B f . 
This implies that cB is a positive linear combination of the rows of B T 
since c f — yA f > 0. Hence cB lies in the polar of {z G M. n ~ d : B T z < 
0} which is the recession cone of P£ proving that the linear program 
min {{—cB)z : z G PJ} is bounded. 

The linear program min {(—cB)z : z G P^} is feasible since is 
a feasible solution. If it is bounded as well then min 

> 0} is feasible and bounded. Hence the dual of 
the latter program max {y ■ b : yA T = c T , yA f < c f } is feasible. This 
shows that a superset of r is a face of A c which implies that r G A c 
since A c is a triangulation. □ 

The reformulations (0) and (Q) imply that G T (b) solves IPA, c {b) if 
and only if both programs have the same optimal solution z* G Z"~ d . 
If G T {b) solves IP A , c (b) then G T (6) also solves IPaM for 

every r 

contained in r. We say that r G A c is associated to IPa,c if for some 
b G NA, G T {b) solves IP A , c ( b ) but ° T \ h ) does not for all' faces r' ^ r 
of A c containing r. Several results about the structure of the subposet 
of faces of A c that are associated to IPa, c can be found in [ 13| . For 



instance, the associated sets of IPa,c occur in saturated chains |13| , 
Theorem 3.1]. For a given b G NA, the most easily solved relaxations 
of IPA, c {b) are those G a (b) where a is a maximal face of A c . We call 
these "top-level" relaxations the Gomory relaxations of IP A,c{b)- 

Definition 2.8. The family of integer programs IPa.c is a Gomory 
family if, for every b G NA, /Pa jC (6) is solved by a group relaxation 
G a (b) where a is a maximal face of the regular triangulation A c . 

Our goal in the rest of this section is to characterize Gomory families 



of integer programs (Theorem 2.17 ). We will assume from now on that 



every integer program in IPa, c has a unique solution which is a stricter 
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notion of genericity of c than requiring A c to be a triangulation p0 |. 
Let O c C N n be the set of all optimal solutions of the programs in 
IPac- The set O c is known to be a down set or order ideal in N n , i.e. 



u E O c and v < u, v G N n implies that u G O c p3[ . For a given 
A, there are only finitely many sets O c as c varies. Two generic cost 
vectors c and c' are equivalent if O c = Q c > and all equivalence classes of 



generic cost vectors are open full dimensional cones in R n [ 2(J . Since c 
is generic, O c is in bijection with NA via the linear map <$a '■ N n — ■> 
where u h-> Am. Let Q u := G M n ~ d : Bz < u, (-cB)z < 0} and 
Qj; := G M n " d : S f z < tt t (u), (-cB)z < 0}. 

Lemma 2.9. (%) ^4 vector u is in O c if and only if Q u fl Z n ~ d = {0}. 
(ii) If u G O c , then the group relaxation G T (Au) solves the integer 
program IPa, c {Au) if and only if Q T U fl Z n ~ d = {0}. 

Proof, (i) The lattice point u belongs to O c if and only if u is the 
optimal solution to IPa, c {Au) which is equivalent to G Z n ~ d being 
the optimal solution to the reformulation (JJ) of IPa, c {Au). Since c 
is generic, the last statement is equivalent to Q u fl Z n ~ d = {0}. The 
second statement follows from the fact that (0) solves (|l|) if and only 
if they have the same optimal solution. □ 



By Lemma |2T9| , it is convenient to use the optimal solution to IPa. c (o) 



as the vector u in (|1]) and (0), and we will do so from now on. For an 
element u G O c and a face r of A c let S(u,r) be the affine semigroup 
u + N(ej : i G r) in N n where denotes the z-th unit vector in R n . 

Lemma 2.10. Suppose u lies in O c . If G T (Au) solves I Pa, c {Au) , then 
G T (Av) solves IPa, c (Av) for all v G S(u,t). 

Proof. If v G S(u,t), then ir T (u) = 7r T (v), and this implies that Q T U = 
Ql. If G T {Au) solves IP A , c (Au), then {0} = Q f u n Z n ~ d = Q; n Z n - d 
for all u G S(u, r). This implies the result by Lemma |2.9| (ii). □ 



Proposition 2.11. The affine semigroup S(u,r) is contained in O c if 
and only ifG T (Au) solves IPa, c {Au). 

Proof. Suppose S(u,t) C O c . Then for all v G S(u,t), Q v = {z G 
Rn -d . ^ < 7r _^j ) B f z < (-cB)z < 0} n Z n - d = {0}. Since 

7Tt=(v) can be any vector in N' T ', and Q T U is bounded by Theorem |2.7| , 
= G : 5 f 2 < 7r r (u), (-cB)z < 0} n Z n ~ d = {0}. Hence, 

by Lemma [T9] (ii), G t (Am) solves JPa )C (Am). 

Conversely, if G T (Au) solves JP Ac (Aw), then {0} = QinZ n ~ d = QIC) 
Z n - d for all v G S(u, r). Since is a relaxation of Q,,, g,,nZ"- d = {0} 
for all v G S(u, r) and hence by Lemma (i), S , (m, r) C C . □ 
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Definition 2.12. For r £ A c and « £ O c , the pair (it, r) is called an 
admissible pair of O c if 

(i) G r (Au) solves IPa, c (Au) or equivalently, S(u,t) C C c , and 

(ii) the support of w is contained in f . 

An admissible pair (u, r) is a standard pair of O c if the affine semigroup 
S(u,t) is not properly contained in another affine semigroup S'(i>,t / ) 
where (w,t / ) is also an admissible pair of O c . 



Definition 2.13. For u £ N n and a face r of A c , we say that the 
polytope Q T U is a standard polytope of JP^ )C if Q 1 u nZ n ~ d = {0} and every 
relaxation of Q T U obtained by removing an inequality in B T z < tt t (u) 
contains a non-zero lattice point. 



Theorem 2.14. The following are equivalent: 

(i) An admissible pair (u,t) is a standard pair of O c 

(ii) The polytope Q T U is a standard polytope of IPa, c - 
(Hi) The face r of A c is associated to IPa, c - 



Proof. The proof of (i) (ii) is the content of Theorem 2.5 in flT3 |. 
The equivalence (i) (Hi) follows from the definition of a standard 
pair and Lemma |2.1U| . □ 

Under the linear map <f>A '■ N n — > NA such that u i— > Au, the affine 
semigroup S(u, r) where (u, r) is a standard pair of O c maps to the 
affine semigroup Au + NA T in NA. Since every integer program in IPa, c 
is solved by one of its group relaxations, O c is covered by its standard 
pairs. We call this cover and its image in NA under cf>A, the standard 
pair decompositions of O c and NA respectively. Since standard pairs 
of O c are determined by the standard polytopes of IPa,c, the standard 
pair decomposition of O c is unique. The terminology used above has 
its origins in |2l| which introduced the standard pair decomposition of a 



monomial ideal. The specialization to integer programming appear in 



19| , §12. D], [0 and ||13|| . For each r £ A c , there are only finitely many 



standard pairs of O c that are indexed by r. Borrowing terminology 



from ||21|| , we call the number of standard pairs of the form (-,r) the 
multiplicity of r in O c . The total number of standard pairs is called the 
arithmetic degree of O c . By Theorem [2.13| , multiplicity of r in O c is the 
number of distinct standard polytopes of IPa,c indexed by f and the 
arithmetic degree of O c is the total number of standard polytopes of 



IPa,c- Lemma [2.15| shows that the maximal faces of A c play a special 
role in the standard pair decomposition of O c . Part (ii) can also be 
deduced from the work of Gomory M. 
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Lemma 2.15. [|T9| , §12. D] (i) (0, a) is a standard pair of O c if and 
only if a is a maximal face of A c . 

(ii) If a is a maximal face of A c , the multiplicity of a in O c is the index 
of the sublattice 7hA a in TLA. 



The index of ZA CT in T,A is also the determinant of A a divided by the 
g.c.d of the maximal minors of A. In contrast to Lemma |2.15| , if r is 
a lower dimensional face of A c , it may not index any standard pair of 
O c . This makes the structure of O c complicated. (See for further 
results.) We give an example below. 



Example 2.16. Consider the following A e Z 3x6 of rank three: 



5 2 1 \ 
5 1 4 2 
5 2 3 / 



The first three columns of A generate cone(A) which is simplicial. If 
c = (21, 6, 1, 0, 0, 0) then the regular triangulation A c is: 



{{1,3,4}, {1,4,5}, {2,5,6}, {3,4,6}, {4,5,6}}. 



The set O c has arithmetic degree 70 which means that O c has 70 stan- 
dard pairs which are listed below. Not all lower dimensional faces of 
A c index standard pairs in this example. 
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standard pairs (-,r) 


{1,3,4} 


(0,-), (es,-), ( e 6,0, (e 5 + e 6 ,-), ( 2e 6, •) 


{1,4,5} 


(0,-), (e2,0> ( e 3,-), ( e 6,-), (e 2 + e 3 ,-), (2e 2 ,-)> 
(3e 2 , •), (2e 2 + e 3 , •) 


{2,5,6} 


(0,-), (e 3 ,-), (2e 3 ,-) 


{3,4,6} 


(0,-), (e 5 ,-), (2e 5 ,-), (3e 5 ,-) 


{4,5,6} 


(0,-), (e 3 ,-), (2e 3 ,-), (3e 3 ,-), (4e 3 ,-) 


{1,4} 


(e 3 + 2e 5 + e 6 , ■), (2e 3 + 2e 5 + e 6 , ■), (2e 3 + 2e 5 , •), 
(2e 3 + 3e 5 ,-), (2e 3 + 4e 5 ,-) 


{1,5} 


(e 2 + e 6 , •), (2e 2 + e 6 , ■), (3e 2 + e 6 , ■) 


{2,5} 


(e 3 + e 4 , •), (64,-), (2e 4 ,-) 


{3,4} 


(e 2 , •), (ei + e 2 , •), (e x + 2e 5 , •), (ei + 2e 5 + e 6 , •), (e 2 + e 5 , •), 


{3,6} 


(e 2 ,-), (e 2 + e 5 ,-) 


{4,5} 


(e 2 + 2e 3 , •), (e 2 + 3e 3 , •), (2e 2 + 2e 3 , •), (3e 2 + e 3 , •), (4e 2j •) 


{5,6} 


(e 2 + 3e 3 ,-) 


{1} 


(e 2 + e 3 + e 6 , ■), (e 2 + e 3 + e 5 + e 6 , •), (e 2 + 2e 6 , ■), 
(e 2 + e 3 + 2e 6 , •), (2e 2 + 2e 6 , •), (e 2 + e 3 + 2e 5 + e 6 , •) 


{3} 


(ei + e 2 + e 6 , ■), (ei + e 2 + 2e 6 , •) 


{4} 


(ei + e 2 + 2e 3 + e 5 , •), (e x + e 2 + 2e 3 + 2e 5 , •), 
(ei + e 2 + 2e 3 + 3e 5 , •), (ei + e 2 + 2e 3 + 4e 5 , ■), 
(ei + 3e 3 + 3e 5 , ■), (e x + 3e 3 + 4e 5 , ■) 


{0} 


(ei + e 2 + 2e 3 + e 5 + e 6 , •), (ei + e 2 + 2e 3 + 2e 5 + e 6 , ■), 
(ei + 2e 2 + e 3 + e 6 , ■), (e x + 2e 2 + e 3 + e 5 + e 6 , ■), 
(ei + 2e 2 + e 3 + 2e 5 + e 6 , •), (ei + 2e 2 + e 3 + 2e 6 , •), 
(ei + 3e 2 + 2e 6 , •) 



Observe that the integer program IP^ c (b) where b = A{e\ + e 2 + e 3 ) 
is solved by G T (b) with r = {1,4,5}. By Proposition [2.3| , Gomory's 
relaxation of IPA, c {b) is indexed by a = {4,5,6} since b lies in the 
interior of the face cone(A a ) of A c . However, neither this relaxation 
nor any nontrivial extended relaxation solves IP^ c (b) since the optimal 
solution ei + e 2 + e 3 is not covered by any standard pair (•, r) where r 
is a non-empty subset of {4,5,6}. 

The results stated thus far give characterizations of Gomory families. 

Theorem 2.17. The following conditions are equivalent: 

(i) IPa, c is a Gomory family. 

(ii) The associated sets of IPa, c are precisely the maximal faces of A c . 
(Hi) (•, r) is a standard pair of O c if and only if r is a maximal face of 
the regular triangulation A c . 

(iv) All standard polytopes of IPa, c are simplices. 
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Proof. The proof follows from Definition [2.8|, Proposition [2.11| and The- 



orem [2.14| . □ 

If there is a generic cost vector c such that for a triangulation A of 
cone(A), A = A c , then we say that A supports the order ideal O c and 
the family of integer programs IPa, c - No regular triangulation of the 
matrix A in Example |2.16| supports a Gomory family. Here is a matrix 



with a Gomory family. 
Example 2.18. Consider the 3x6 matrix 

/ 1 1 1 1 1 
A= 1 1 1 2 2 
\ 1 2 3 4 

In this case, cone(A) has 14 distinct regular triangulations and 48 dis- 
tinct sets O c as c varies among all generic cost vectors. Ten of these 
triangulations support Gomory families; one for each triangulation. For 
instance, if c = (0, 0, 1, 1, 0, 3), then 

A c = {a 1 = {1,2,5}, a- 2 = {1,4, 5}, cr 3 = {2,5,6}, a 4 = {4,5,6}} 

and IPa,c is a Gomory family since the standard pairs of O c are: 
(0,Oi), (e 3 ,<7i), (e 4 ,cri), (0,a 2 ), (0,er 3 ), and(0,<r 4 ). 

3. Total dual integrality and Gomory families 

We now relate the notion of total dual integrality [fT6"| , §22] to Gomory 
families. Recall that TLA = Z d by assumption. 

Definition 3.1. The system yA < c is totally dual integral (TDI) if 
LPa,c(P) has an integral optimal solution for each b G cone(A) n Z° '. 

Definition 3.2. The regular triangulation A c is unimodular if 1A a = 
Z d for every maximal face a G A c . 

Theorem 3.3. T7ie system yA < c is TDI if and only if the regular 
triangulation A c is unimodular. 

Proof. The regular triangulation A c is the normal fan of P c by Propo- 
sition |2.4j , and it is unimodular if and only if ZA a = 7L d for ev- 
ery maximal face o G A c . This is equivalent to saying that every 
o G cone (A a ) n Z d lies in NA a for every maximal face a of A c . By 
Lemma |2.3| , this happens if and only if LP^ c (b) has an integral opti- 
mum for all b G cone(A) fl 7L d . □ 

Corollary 3.4. If yA < c is TDI then IPac is a Gomory family. 
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Proof. By Lemma |2.15| , (0, o~) is a standard pair of O c for every maximal 



face a of A c . Theorem |3]^ implies that coroe(A) is unimodular (i.e. 
ZA a = Z d ), and therefore NA a = cone(Acr) flZ d for every maximal face 
cr of A c . Hence the semigroups NA a arising from the standard pairs 
(0, a) as a varies over the maximal faces of A c cover NA Therefore 
the only standard pairs of O c are (0, a) as o varies over the maximal 
faces of A c . The result then follows from Theorem [2. 17] (ii). □ 



When yA < c is TDI, the multiplicity of any maximal face a of A c 
in O c is one, and all other faces have multiplicity zero. While this is 
sufficient for IPa, c {o) to be a Gomory family, it is far from necessary. 
TDI-ness guarantees that LPa, c {o) has an integral optimum for every 
integral b in cone(A). In contrast, if IP^ c (b) is a Gomory family, the 
linear optima of the programs in LPa, c may not be integral. 

If A is unimodular (i.e., ZA = ^ d for every nonsingular maximal 
submatrix A of A), then the feasible regions of the linear programs 
in LPa, c have integral vertices for each integral b G cone(A) R Z d , 
and yA < c is TDI for all c. Hence if A is unimodular, then IPa, c 
is a Gomory family for all generic cost vectors c. However, just as 
integrality of the optimal solutions of programs in LPa, c is not necessary 
for IPa, c to be a Gomory family, unimodularity of A is not necessary 
for IPa, c to be a Gomory family for all c. 

Example 3.5. Consider the seven by twelve integer matrix 

/lOOOOOlllllO^ 
010000110001 
001000101001 
A = 1 1 1 
000010001010 
000001000111 

V000000111111/ 

of rank seven. The maximal minors of A have absolute values zero, 
one and two and hence A is not unimodular. This matrix has 376 
distinct regular triangulations supporting 418 distinct order ideals O c . 
In each case, the standard pairs of O c are indexed by just the maximal 
simplices of the regular triangulation A c that supports it. Hence IPa, c 
is a Gomory family for all generic c. 

4. A-NORMAL MATRICES 

In Section 3 we saw that unimodularity of A or more generally, uni- 
modularity of a regular triangulation of cone(A), gives rise to Gomory 
families of integer programs. In this section, we identify a larger set of 
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matrices and cost vectors that give rise to Gomory families. A com- 
mon property of unimodular matrices and matrices with a unimodular 
triangulation is that they form a Hilbert basis for cone(A). In other 
words, NA equals cone(A) D Z d for such matrices. Borrowing a term 
from commutative algebra we make the following definition. 

Definition 4.1. A d x n integer matrix A is normal if the semigroup 
NA equals cone(A) n Z d . 

We first note that if A is not normal, then IPa, c need not be a 
Gomory family for any cost vector c. 

Example 4.2. The non-normal matrix A = \ \ } \ \\ gives rise 



x 1 3 4 / 

to 10 distinct order ideals O c supported on its four regular triangula- 
tions{{l,4}},{{l,2},{2,4}},{{l,3},{3,4}}and{{l,2},{2,3},{3,4}}. 
Each O c has at least one standard pair that is indexed by a lower di- 
mensional face of A r . 



The matrix in Example |2.16| is also not normal and has no Gomory 



families. These examples show that normality of A is necessary for the 
existence of Gomory families. However, we do not know at this time 
whether every normal matrix A has some generic cost vector c such 
that IPa.c is a Gomory family. Our goal is to show that under certain 
additional conditions, normal matrices do give rise to Gomory families. 

Definition 4.3. A d x n integer matrix A is A-normal if it has a 
triangulation A such that for every maximal face a G A, the columns 
of A in cone(A a ) form a Hilbert basis for cone(A a ). 

Remark 4.4. If A is A-normal for some triangulation A, then it is 
normal. To see this note that every lattice point in cone(A) lies in 
cone(A a ) for some maximal face a G A. Since A is A-normal, this 
lattice point also lies in the semigroup generated by the columns of A 
in cone(Aj) and hence in NA 

Observe that A is A-normal with respect to all the unimodular tri- 
angulations of cone(A). Hence triangulations A with respect to which 
A is A-normal generalize unimodular triangulations of cone(A). 



Examples |4.5| and [4.6| show that the set of matrices where cone(A) 



has a unimodular triangulation is a proper subset of the set of A-normal 
matrices which in turn is a proper subset of the set of normal matrices. 

Example 4.5. Examples of normal matrices with no unimodular tri- 
angulations can be found in [0 and |7|]. If cone(A) is simplicial for 
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such a matrix, A will be A-normal with respect to its coarsest (regu- 
lar) triangulation A consisting of the single maximal face with support 
cone(A). For instance, consider the following example taken from [|7|]: 



Here cone(A) has 77 regular triangulations and no unimodular trian- 
gulations. Since cone(A) is simplicial, A is A-normal with respect to 
its coarsest regular triangulation {{1,2,3,8}}. 

Example 4.6. There are normal matrices A that are not A-normal 
with respect to any triangulation of cone(A). To see such an example, 
consider the following modification of the matrix in Example [4.5| that 
appears in |n§ Example 13.17] : 



This matrix is again normal and each of its nine columns generate an 
extreme ray of cone(A). Hence the only way for this matrix to be 
A-normal for some A would be if A is a unimodular triangulation of 
cone(A). However, this cone(A) has no unimodular triangulations. 

Theorem 4.7. If A is A-normal for some regular triangulation A then 
there exists a generic cost vector c G U 1 such that A = A c and IPa, c 
is a Gomory family. 

Proof. Without loss of generality we can assume that the columns of A 
in cone(A a ) form a minimal Hilbert basis of this cone for any maximal 
face a of A. If there were a redundant element, the smaller matrix 
obtained by removing this column from A would still be A-normal. 

For a maximal face a G A, let o~i n C {1, . . . , n} be the set of indices of 
all columns of A lying in cone{A a ) that are different from the columns 
of A a . Suppose a^,... ,a ik are the columns of A that generate the 
one dimensional faces of A, and d G M™ a cost vector such that A = 
A c /. We modify d to obtain a new cost vector c G M™ such that 
A = A c as follows. For j = 1, . . . , k, let c ij := d i . If j G a in for 
some maximal face a G A, then aj = ^2 iea < \ < 1 and 
we define Cj := ^j go .AjCj. Hence, for all j G a in , (aj,Cj) G R d+1 lies 



.4 



/ 1 1 1 1 1 1 \ 
10 112 2 2 
1 1 2 2 3 3 

^00012345/ 



.4 



/ 1 1 1 1 1 1 \ 
1 1 1 2 2 2 
1 1 2 2 3 3 
000012345 

y 1 1 1 1 1 1 1 1 1 J 
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in C a := cone((ai, Ci) : % G o) — cone((a,i, d : i G cr) which was a 
facet of C = cone((ai, c'j) : z = 1, . . . , n). If y G M d is a vector as in 



Definition |2.1| showing that cr is a maximal face of A c > then y • ai = Ci for 
alH G crU(7j n and y-cij < Cj otherwise. Since cone(A-) = cone(A aUlJin ) , 
we conclude that cone(A a ) is a maximal face of A c . 

If b G NA lies in cone(A a ) for a maximal face a G A c , then IP^ c (b) 
has at least one feasible solution u with support in cr U crj n since A 
is A-normal. Further, (b, c ■ u) = (Au, c ■ u) lies in C a and all feasible 
solutions of IPa,c(P) with support in aUa. in have the same cost value by 
construction. Suppose v G N n is any feasible solution of IPA,c(b) with 
support not in a U a in . Then c - u < c - v since (ctj, ci) G C CT if and only 
if i G cr U a in and is a facet in the /ou>er envelope of C. Hence the 
optimal solutions of IPa, c (P) are precisely those feasible solutions with 
support in cr U a in . The vector b can be expressed as b = b' + X^e^ z i a i 
where G N are unique and b' G {J2iea -^ a » : < Aj < 1} fl Z d is also 
unique. The vector b' = J2jea in r j a j wri ere rj G N. Setting Ui = Zi for 
all i G cr, itj = r_j for all j G crj„ and li^ = otherwise, we obtain all 
feasible solutions u of IPA, c (b) with support in cr U a in . 

If there is more than one such feasible solution, then c is not generic. 
In this case, we can perturb c to a generic cost vector c" = c + euj 
by choosing 1 3> e > 0, Uj <^ whenever j = i\, . . . ,%k and Uj = 
otherwise. Suppose U\, . . . ,u t are the optimal solutions of the integer 
programs IP A ,c"(b') where V G {E ie<J : < Aj < 1} fl Z d . (Note 
that t = |{Ei G( r^ a i : < Ai < 1} n Z d | is the index of ZA a in 
ZA) The support of each such ii« is contained in o~i n . For any 6 G 
cone (Aj) H Z d , the optimal solution of IPA, c "{b) is hence u = Ui + z 
for some i G {1, . . . , i} and z G N n with support in a. This shows 
that KL4 is covered by the affine semigroups <pA(S(ui, a)) where a is a 
maximal face of A and as above for each a. By construction, the 
corresponding admissible pairs (itj,cr) are all standard for O c ii. Since 
all data is integral, c" G Q n and hence can be scaled to lie in Z n . 
Renaming c" as c, we conclude that IPa, c is a Gomory family. □ 

Corollary 4.8. Let A be any normal matrix such that cone(A) is sim- 
plicial, and let A be the coarsest triangulation whose single maximal 
face has support cone(A). Then there exists a cost vector c G Z n such 
that A = A c and IPa, c is a Gomory family. 



Example 4.9. Consider the normal matrix in Example [2.18| . Here 



cone(A) is generated by the first, second and sixth columns of A 
and hence A is A-normal with respect to the regular triangulation 
{{1, 2, 6}}. There are 13 distinct sets O c supported on A. Among the 
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13 corresponding families of integer programs, only one is a Gomory 
family. A representative cost vector for this IPa,c is c = (0, 0, 4, 4, 1, 0). 
The standard pair decomposition of O c is the one constructed in The- 



orem [4.7| . The affine semigroups S(-,a) from this decomposition are: 

S(0,a), S(e 3 ,a), S , (e 4 ,cr), and S(e 5 , a). 

Note that A is not A-normal with respect to the regular triangulation 
supporting the Gomory family IPa, c in Example p. 18 , The columns of 



A in cone(A ai ) are the columns of A ai and A 3 . The vector (1, 2, 2) is in 
the minimal Hilbert basis of cone(A ai ) but is not a column of A. This 
example shows that a regular triangulation A of cone(A) can support 
a Gomory family even if A is not A-normal. The Gomory families in 
Theorem |4.7| have a very special standard pair decomposition. 



5. Hilbert covers and Gomory families 

The results in the previous section lead to the following problem. 

Problem 5.1. If A G Z dxn is a normal matrix, does there exist a 
generic cost vector c G Z n such that IPa,c is a Gomory family? 

We do not know the answer to this question. However, in this section 
we answer a stronger version of this question for small values of d and 
state our observations for general d. We begin with the following result. 

Theorem 5.2. If A G Z dxn is a normal matrix and d < 3, then there 
exists a generic cost vector c G Z n such that IPa,c is °> Gomory family. 

Proof. It is known that if d < 3 then cone(A) has a regular unimodular 
triangulation A c |T7[. The result then follows from Corollary |3~3. □ 



Before we proceed, we rephrase Problem |5J] in terms of covering 
properties of cone(A) and NA along the lines of 0, ||, M, |7J and 
171 . To obtain the same set up as in these papers we assume in this 



section that A is normal and the columns of A form the unique minimal 
Hilbert basis of cone(A). Using the terminology in ||, the free Hilbert 
cover problem asks whether there exists a covering of NA by semigroups 
NA T where the columns of A T are linearly independent. The unimodu- 
lar Hilbert cover problem asks whether cone(A) can be covered by full 
dimensional unimodular subcones cone(A T ) (i.e., ZA T = Z d ), while the 
stronger unimodular Hilbert partition problem asks whether cone(A) 
has a unimodular triangulation. (Note that if cone(A) has a unimod- 
ular Hilbert cover or partition using subcones cone(A T ), then NA is 
covered by the semigroups NA,-.) All these problems have positive an- 
swers if d < 3 since cone(A) admits a unimodular Hilbert partition in 
this case [0], JTTj. Normal matrices (with d = 4) such that cone(A) has 
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no unimodular Hilbert partition can be found in and [J/J. Examples 
(with d — 6) that admit no free Hilbert cover and hence no unimodular 
Hilbert cover can be found in || and |4j]. 

When < c is TDI, the standard pair decomposition of NA induced 



by c gives a unimodular Hilbert partition of cone(A) by Theorem pT3. 



An important difference between Problem ETT] and the Hilbert cover 



problems is that affine semigroups cannot be used in Hilbert covers. 
Moreover, affine semigroups that are allowed in standard pair decom- 
positions come from integer programming. If there are no restrictions 
on the affine semigroups that can be used in a cover, NA can always 
be covered by full dimensional affine semigroups: for any triangulation 
A of cone(A) with maximal subcones cone(A a ), the affine semigroups 
b + NA a cover NA as b varies in ^ '■ < A; < 1} H Z d and 

cr varies among the maximal faces of the triangulation. A partition of 
NA derived from this idea can be found in |18|, Theorem 5.2]. 

In order to state our main theorem, we recall the notion of super- 
normality which was introduced in Jll| . 



Definition 5.3. A matrix A e % dxn j s supernormal if for every sub- 
matrix A' of A, the columns of A that lie in cone(A') form a Hilbert 
basis for cone(A'). 

Proposition 5.4. For A 6 Z dxn , the following are equivalent: 

(i) A is supernormal, 

(ii) A is A-normal for every regular triangulation A of cone(A), 

(iii) Every triangulation of cone(A) in which all columns of A generate 
one dimensional faces is unimodular. 



Proof. The equivalence of (i) and (iii) was established in O, Propo- 



sition 3.1]. Definition shows that (i) =>• (ii). Hence we just need 
to show that (ii) =>■ (i). Suppose that A is A-normal for every regular 
triangulation of cone(A). In order to show that A is supernormal we 
only need to check submatrices A' where the dimension of cone(A') is 
d. Choose a cost vector c with q ^> if the i-th column of A does not 
generate an extreme ray of cone(A'), and q = otherwise. This gives 
a polyhedral subdivision of cone(A) in which cone(A') is a maximal 
face. There are standard procedures that will refine this subdivision 
to a regular triangulation A of cone(A). Let T be the set of maximal 
faces a of A such that cone(A a ) lies in cone(A'). Since A is A-normal, 
the columns of A that lie in cone{A a ) form a Hilbert basis for cone(A a ) 
for each a G T. However, since their union is the set of columns of A 
that lie in cone(A'), this union forms a Hilbert basis for cone(A'). □ 
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It is easy to catalog all A-normal and supernormal matrices, of the 
type considered in this paper, for small values of d. We say that the 
matrix A is graded if its columns span an affine hyperplane in M. d . If 
d — 1, cone(A) has n triangulations {{i}} each of which has the unique 
maximal subcone cone(Ai) whose support is cone(A). If we assume that 
a i < a 2 < ' ' • < o-ni then A is normal if and only if either a\ = 1, or 
a n = —1. Also, A is normal if and only if it is supernormal. If d = 2 
and the columns of A are ordered counterclockwise around the origin, 
then A is normal if and only if det(ai, aj + i) = 1 for all % = 1, . . . ,n — l. 
Such an A is supernormal since it is A-normal for every triangulation 
A — the Hilbert basis of a maximal subcone of A is precisely the set 
of columns of A in that subcone. If d = 3 then as mentioned before, 
cone(A) has a unimodular triangulation with respect to which A is A- 
normal. However, not every such A needs to be supernormal: we saw 
that the matrix in Example [2.1S| is not A-normal for the A supporting 
the Gomory family in that example. If d — 3 and A is graded, then 
without loss of generality we can assume that the columns of A span the 
hyperplane x\ = 1. If A is normal as well, then its columns are precisely 
all the lattice points in the convex hull of A. Conversely, every graded 
normal A with d = 3 arises this way — its columns are all the lattice 
points in a polygon in M? with integer vertices. In particular, every 
triangulation of cone(A) that uses all the columns of A is unimodular. 
Hence, by Proposition |5.4j , A is supernormal, and therefore A-normal 
for any triangulation of A. 

Theorem 5.5. Let A e Z dxn be a normal matrix of rank d. 

(i) If d = 1,2 or A is graded and d = 3, every regular triangulation 
of cone(A) supports at least one Gomory family. 

(ii) If d = 2 and A is graded, every regular triangulation of cone(A) 
supports exactly one Gomory family. 

(iii) If d = 3 and A is not graded, or if d = 4 and A is graded, then 
not all regular triangulations of cone(A) may support a Gomory 
family. In particular, A may not be A-normal with respect to 
every regular triangulation. 

Proof, (i) If d — 1, 2 or A is graded and d — 3, A is supernormal and 



hence by Proposition |5.4j and Theorem |4.7| , every regular triangulation 
of cone(A) supports at least one Gomory family. 

(ii) If d — 2 and A is graded, then we may assume that 



.4 
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In this case, A is supernormal and hence every regular triangulation 
A of cone(A) supports a Gomory family by Theorem [4.7| . Suppose 
the maximal cones of A, in counter-clockwise order, are Ci,... ,C r . 
Assume the columns of A are labeled such that Cj = cone(cii-i, a %) f° r 
i — 1, . . . ,r, and the columns of A in the interior of Cj are labeled in 
counter-clockwise order as bn, . . . , b^.. Hence the n columns of A from 
left to right are: 

a o, bn, ■ ■ ■ , fri/tp ai, &21) • • • ; &r-ij b r i, ■ ■ • , b r k r , a r - 

Indexing the columns of A by their labels, the maximal faces of A are 
(jj = {i — 1, i} for i = l,... ,r. Let be the unit vector of M™ indexed 
by the true column index of a, in A and be the unit vector of M" 
indexed by the true column index of b{j in A. Since the columns of 
A form a minimal Hilbert basis of cone(A), is the unique solution 
to JP J 4 )C (aj) for all c and eij is the unique solution to 7P J 4, C (%) for all 



c. Hence the standard pairs of Theorem fO] are (0, cij) and (e^, <7j) for 
i = 1, ... ,r and j = 1, . . . , 

Suppose A supports a second Gomory family IPa,u- Then every 
standard pair of 0^ is also of the form (•, a"j) for a,i G A, and r of them 
are (0,(Tj) for i = 1, . . . , r. The remaining standard pairs are of the 
form (e^, cr^). To see this, consider the semigroups in KL4 arising from 
the standard pairs of O w . The total number of standard pairs of O c 
and O w are the same. Since the columns of A all lie on x\ = 1, no two 
bijS can be covered by a semigroup coming from the same standard pair 
and none of them are covered by a semigroup (0, <7j). We show that if 
(e^, is a standard pair of O w then k = i and thus O w = O c . 

If r = 1, the standard pairs of O w are (0, a"i), (en, a\), . . . , (e^, o"i) 
as in Theorem fO. If r > 1, consider the last cone C r = cone(a r ^\, a r ). 



If a r _i is the second to last column of A, then C r is unimodular and 
the semigroup from (0, a r ) covers C r fl Z 2 . The subcomplex comprised 
of C\ , . . . , C T — i is a regular triangulation A' of cone(A') where A' is 
obtained by dropping the last column of A. Since A' is a normal graded 
matrix with d = 2 and A' has less than r maximal cones, the standard 
pairs supported on A' are as in Theorem [4.7| by induction. If a T _i is 
not the second to last column of A then b r ^ r , the second to last column 
of A is in the Hilbert basis of C T but is not a generator of C r . So O w 
has a standard pair of the form (e r k r , (?i)- If o-i ^ o r , then the lattice 
point b r k r + a r cannot be covered by the semigroup from this or any 
other standard pair of O w . Hence Oi = a r . By a similar argument, the 
remaining standard pairs indexed by a T are (e r ^ r -i), cr r ), . . . , (e r ±, a r ) 
along with (0, a r ). These are precisely the standard pairs of O c indexed 
by a r . Again we are reduced to considering the subcomplex comprised 
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of Ci, . . . , C r _i and by induction, the remaining standard pairs of O v 
are as in Theorem 



(iii) The 3x6 normal matrix A of Example |2.18| has 10 distinct Gomory 
families supported on 10 out of the 14 regular triangulations of cone(A). 
Furthermore, the normal matrix 



A 



( 1 
1 




i 
i 

2 
4 



1 \ 






has 11 distinct Gomory families supported on 11 out of the 19 regular 
triangulations of cone(A). □ 



6. Computations and Commutative Algebra 

As mentioned in the introduction, the computations in this paper 
were done using the connections of this material to commutative alge- 



bra (see [191], pll and |23[). In this section we give a brief description 



of our methods. 

Our general feeling is that Problem 5T has a negative answer. To 
check whether a matrix A has a Gomory family, we first need to com- 
pute all the distinct sets of optimal solutions O c (to the programs in 
IPa,c) that arise as c varies among the generic cost vectors with re- 
spect to A. As mentioned in Section 2, there are only finitely many 
such sets for a fixed A. To check whether IPa, c is a Gomory family, 
we need to compute the standard pair decomposition of O c : IPa, c is a 
Gomory family if and only if all the standard pairs of O c are indexed 
by maximal faces of A c . 

A monomial x u in the polynomial ring S := Q[xi, . . . ,x n ] is the 
product x u = x" 1 • • -x^ n where u = (ui, . . . ,u n ) G N n . An ideal M 
in S is called a monomial ideal if it is generated by monomials. Since 
every ideal in S is finitely generated, M = (x vi , ... ,x Vt ) for a set of 
minimal generators X j • • • i X • The toric ideal of A, denoted as I a is 
the binomial ideal in S defined as: 

(x u - x v : u, v G N n and Au = Av). 

The cost of a monomial x u with respect to a cost vector c G M. n is the 
dot product c-u and the initial term of a polynomial f = ^u% u G S is 
the sum of all terms in / of highest cost. For any ideal I G S, the initial 
ideal of I with respect to c, denoted as in c (I), is the ideal generated 
by all the initial terms of all polynomials in I. These concepts come 
from the theory of Grobner bases for polynomial ideals [H]. The toric 
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ideal I a provides the algebraic link between integer programming and 
Grobner basis theory. For an introduction to this connection see |L9| . 

Proposition 6.1. (%) A cost vector c G M. n is generic for A if and only 
if the initial ideal inJjA) is a monomial ideal. 

(ii) For a generic c, the vector u belongs to O c if and only if x u is not 
in the intial ideal in c (lA)- O 

There are only finitely many distinct initial ideals for a polynomial 



ideal fL9| , and hence there are only finitely many distinct sets O c as c 
varies among the generic cost vectors. We say that two cost vectors c 
and d in R n are equivalent if in c (l a) = in c '(lA)- 

Theorem 6.2. [ pC| , Theorem 3.10] (i) Each equivalence class of generic 
cost vectors form an open full dimensional polyhedral cone in W 1 . The 
closure of this cone is called the Grobner cone of c where c is any vector 
in the interior of the cone. 

(ii) The collection of all Grobner cones of A form a complete polyhedral 
fan in R™ called the Grobner fan of A. 

(Hi) The Grobner fan of A is the normal fan of an (n — d)- dimensional 
polytope called the state polytope of A. □ 

By the above results, O c can be computed implicitly by computing 
the monomial ideal in c (lA)- This can be done using a computer algebra 
package like Macaulay 2 [IIJ. In order to find all initial ideals of I a, 
we use the software package TiGERS for enumerating the vertices 
of the state polytope of A. At each vertex, TiGERS returns the initial 
ideal induced by a vector in the interior of the normal cone at that ver- 
tex. The standard pair decomposition of the set of monomials outside 
a monomial ideal described in terms of its minimal generators can be 



calculated using Macaulay 2. See the chapter Monomial Ideals in |12|. 

To obtain a normal matrix of the type discussed in this paper, it suf- 
fices to start with an arbitrary set of vectors a 1; . . . , a p G Z d such that 
cone(ai, . . . , a p ) is pointed and full dimensional and then to compute 
the Hilbert basis of cone(ai, . . . , a p ). The elements in the Hilbert basis 
form the columns of a normal matrix. We used the package Normaliz 
by Bruns and Koch || to compute Hilbert bases. 

The regular triangulation A c of cone(A) is a pure d-dimensional com- 
plex of cones and hence every maximal face a G A c has cardinality d. 
Hence we get the following algorithm. 

Algorithm 6.3. How to check if A G Z dx?l gives a Gomory family. 

(i) Compute all the initial ideals of the toric ideal I a using TiGERS. 

(ii) For each initial ideal in c {lA)'- 

(a) Use Macaulay 2 to find its standard pairs. 
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(b) IPa,c is a Gomory family if and only if every set r in a standard 
pair (-,r) has cardinality d. 

A monomial x u is square-free if u G {0, l} n and a monomial ideal 
M is square-free if all its minimal generators are square-free. The 
radical of a monomial ideal M in 5 is the ideal vM := (/ : / r G 
M for some r G N). The radical yM is a square-free monomial ideal. 
The Stanley Reisner ideal of the regular triangulation A c is the square- 
free monomial ideal 

(Hi<z T Xi : r is a minimal non-face of A c ) . 



Theorem 6.4. [ |19| , Theorem 8.3] The Stanley- Reisner ideal of the reg- 
ular triangulation A c is the radical ofin c (lA)- □ 

If two distinct initial ideals in c {lA) and in^^A) have the same rad- 
ical, then A c = A c > and we say that A c supports these initial ideals. 
Several initial ideals of I a may have the same radical. Using the above 
theorem, the initial ideals of I a output by TiGERS can be grouped 
according to their radicals or equivalently, the regular triangulations 
supporting them. This, in combination with Algorithm |S.3|, allows us 
to check whether a regular triangulation supports a Gomory family. 



Recall from Theorem |3.3| that yA < c is TDI if and only if A c is 
unimodular. Using the following result, we obtain an algebraic check 
for TDI-ness of yA < c. 

Theorem 6.5. Corollary 8.9] The regular triangulation A c is uni- 
modular if and only if the initial ideal in c (lA) is square-free. □ 
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